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We study an integrable spin chain with an alternating array of spins S = 1/2, 1 in external 
magnetic fields using the Bethe ansatz exact solution. The calculated magnetization possesses a 
cusp structure at a critical magnetic field H = He , at which the specific heat shows a divergence 
property. We also calculate finite-size corrections to the energy spectrum, and obtain the critical 
exponents of correlation functions with the use of conformal field theory (CFT). Low-energy 
properties of the model are described by two c — 1 U(l) CFTs in H < He and one c = 1 17(1) 
CFT in H > H c . 
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It is now well established that the antiferromagnetic 
quantum spin chains with half-integer spins have mass- 
less excitations, whereas those with integer spins have 
massive ones (Haldane gap), although the ground state 
of both cases are known to be singlet. Experimentally, 
another interesting class of quantum spin systems have 
been found, i.e. the spin systems with an alternating 
array of integer and half-integer spins both for ferromag- 
netic and antiferromagnetic casear. The systems found 
so far have the magnetically ordered ground state even 
for the antiferromagnetic case, for which quantum fluc- 
tuations may not play a crucial role. It is interesting to 
explore what happens if the magnetic order disappears, 
for example, in the presence of certain interactions be- 
yond nearest neighbor ones. If such a system is realized, 
quantum fluctuations are expected to become important. 
This problem may thus provide us with a new interesting 
paradigm of quantum spin chains. 

As a first step to address this problem of alternating 
spin chains with a singlet ground state, itjnA¥.i)e inter- 
esting to study an integrable spin modelB Beth' whose 
ground state is singlet. In this paper, we study the inte- 
grable alternating spin model of S = 1/2, 1, with a spe- 
cial attention to its critical properties in magnetic fields. 
We systematically evaluate the magnetization, the spe- 
cific heat and the critical exponents of correlation func- 
tions, and discuss their characteristic properties. 

We consider a system of N spins • • 

of spin-i and N spins Si, S3, ■ • • , S2N-1 of spin-1 in an 
externaljpagnetic field H, the Hamiltonian for which is 
given byo ) 
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nian contains rather complicated next-nearest-neighbor 
interactions in addition to nearest-neighbor Heisenberg- 
type interactions, in order for the model to be inte- 
grable^. TheL ground state of the model is known to be 
singletoHotP , for which, as will be shown below, the 
antiferromagnetic correlation is dominant. In spite of 
its ugly appearance (1), we wish to note that this sys- 
tem may provide a starting model to discuss fundamen- 
tal properties of alternating spin chains with quantum 
disordered (singlet) ground states. In the end of the pa- 
per, we shall give some discussions about what would be 
expected beyond this special model. 

Let us first summarize the basic equations necessary 
for the following analysis. The Bethe equation obtained 
ae quantum inverse scattering method for (1) is 
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(2) 



where L = 2N is assumed to be a multiple of 4. Here 
we have taken the fully polarized state as a reference 
state, and denotes the number of "particles" by M, which 
are created by the spin-lowering operator. So, the z- 
component of the total spin is to be S z — 3N/2 — M. 
One can see two kinds of phase factors reflecting the 
mixture of spin 1 and 1/2 in the left-hand side of (2). 
The total energy is 
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E = -N 
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where periodic boundary conditions are imposed: cr 2r j = 
C2n+2JV an< ^ Szn—i = S2 n -i+2N ■ Note that the Hamilto- 



and the total momentum is given by P — Im Y)j^~ 1 log[(Aj - 
f)(Aj — *)]• So far, thermodynamic properties of the 
model have been investigated, revealing some properties 
thi6£6(&f) the ordinary Heisenberg spin chain 
Also, conformal properties have been discussed 
to some extent: the system at zero field is classified as 
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CFT of central charge c = 20' p 1 , and is conjectured to 
be-described by level-1 SU(2) Wess-Zumino-Witten mod- 
elaf. In what follows, we systematically study the model 
in the presence of magnetic fields, and show that below 
a certain critical field He, low-energy properties of the 
model are described by two c = 1 U(l) CFTs, while for 
H > H c by one c = 1 U(l) CFT. Accordingly, the sys- 
tem shows quite different critical behaviors for low fields 
H < H c and high fields H > H C - 

We start with the calculation of the magnetization 
curve. To this end we need the Bethe equations in the 
thermodynamic limit. For large enough L, the ground 
state properties can be described by Mr real rapidities 
A R and 2-string complex rapidities A", 



Xf + -i(3-2a) (a = 1,2) 



(4) 



where A^p (j = 1, 2, • • • , Mc) rsB xeal number which la- 
bels the center of the 2-stringatJacJ. The Bethe equa- 
tion (2) is thus converted into integral equations for the 
density functions p R (A R ) and p (A ) for these rapidi- 
ties 
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where the integral kernels are given by -Krr(A) = K\{X), 
K CC (X) = {2K 1 (X)+K 2 {X)}, K CR (X) = K KC {\) = 
{K 1 / 2 (X) + K 3 / 2 (X)}, respectively, in terms of the func- 
tion K n {x) = (n/Tr){n 2 + x 2 }~ 1 . The driving terms read 



A, 



K x + K 1/2 
K x + K 3 



and the symbol ® denotes the convolution integral in 
the finite range of [A~,A+] respectively for the sectors 
of a=R, C. 

By-jpalculating the above integral equations (5) numer- 
icallytP, we have computed the magnetization A4(= S z ) 
as a function of the magnetic field. The results are shown 
in Fig.l. We can see that the magnetization begins lin- 
early in magnetic fields (finite spin susceptibility) since 
the ground state is singlet, and increases smoothly up 
to the critical field Ha., where we encounter a charac- 
teristic cusp structureLP. The existence of the cusp at 
H = He indicates that the magnetization for H < He 
consists of two contributions, i.e. M. = A^r + A^c and 
one of them, A^c, saturates at H = He- For H > He, 
the remaining massless spin excitations still contribute 
to the magnetization, and then at H = 6 the system 
is completely polarized. These characteristic properties 
directly reflect that the system consists of two kinds of 
different spins. The anomaly at H — He also appears in 
the specific heat. We show the computed results for the 
coefficient 7 of the T-linear specific heat in Fig. 2. The 
divergence property for H — > He implies that the ve- 
locity for one of spin excitations tends to zero when the 
system approaches the critical field He , since the specific 
heat coefficient 7 is inversely proportional to the velocity. 
To see this point clearly, we have plotted two velocities 



wr and vq in Fig. 3. It is seen that vq indeed vanishes 
at H = He- Therefore, we can say that there exist 
two species of spin excitations which behave differently 
in magnetic fields. They have the distinct pseudo-Fermi 
surfaces pp and pp, both of them decrease with the in- 
crease of magnetic fields. In particular, at H = He, Pp 
vanishes and the corresponding spin excitation becomes 
massive. In what follows we refer to these spin excita- 
tions as spinons in the R- and C-sectors respectively. 

In order to further study critical properties of the 
model in magnetic fields, we now calculate the finite-size 
spectrum. Following a standard techniquetP we obtain 
the corrections to the energy 
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L v " L 
where £oo is the bulk energy density in the thermody- 
namic limit. Also, the finite-size corrections to the mo- 
mentum are 



P= V [-2p^AD a - T (A+-A-)} (8) 



a=R,C 



where "the Fermi momenta" for two spinons are given 
by pp = 7r/4 - ttMr/L and pp = 7r/2 - nMc/L with 
A4 R = L/4 - Mr and M c = L/2 - 2M C . Here AD a 
labels an excitation across the two Fermi points, which 
carries the large momentum transfer 2ppAD a . 

We note that both of the above expressions (7) and 
(8) completely fit in with the finite-size scaling formu- 
lae in CFT@). Namely, the corrections —nv a /(6L) to 
the ground state energy in (7) indicate that two kinds of 
spinons are described by two c = 1 CFTs. Also, the con- 
formal dimensions A^ which enter the excitation energy 
as well as the momentum are of desirable forms predicted 
by CFT, 
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forj2L= R, C. Here the so-called dressed charge matrix 
tiajEl 1 is determined by the following integral equations, 



(10) 



where K is the 2x2 matrix of integral kernels which 
has appeared in eq.(5), and / is the identity matrix. 
The above dressed charge matrix determines the criti- 
cal exponents of the present model in magnetic fields. 
Note that the conformal dimensions (9) both for R- and 
C-sectors have a typical form inherent in c = 1 Gaus- 
sian CFT with U(l) symmetry, which has previously ap- 
peared in various contexts, for instance, in the analy- 
sis of interacting electron systemsliiP. Recall here that 
for H > He, spinons of the C-sector become massive 
with vanishing vc ■ Therefore we can see that low-energy 
properties of the present model in magnetic fields are de- 
scribed by two c = 1 CFTs for < H < He and one 
c = 1 CFT for H > H c . 

Having specified the critical behavior in magnetic 
fields by CFT, we now read critical exponents of cor- 
relation functions from the conformal dimensions (9) by 
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suitably choosing the quantum numbers. For this pur- 
pose, we should note that the quantum numbers in (9) 
are subject to the selection rules 



spin-spin relaxation rates 1/T 2 g from eq. (13)0) 



AAf R 



ADr 



AMr 



(mod 1) , (11) 



2 ' 2 
which directly reflect that spinons can be regarded as 
interacting bosons. 

For H < He both of two kinds of spinons control the 
long-distance behavior of correlation functions, since vr 
and vc have finite values. Thus choosing the quantum 
numbers as 



(AMr.AMc.ADr.ADc) 



(0,0,1,0) 
(0,0,0,1) 



(12) 



we find the spin-spin correlation function to have the 
asymptotic form, 

ao a 1 cos{2p^:x) a 2 cos(2p%x) 
Xz{x) ~ — + -r + ^5 (13) 

X z X a ° X a s 

Note that the spin-spin correlation function has two os- 
cillating terms, reflecting the presence of two-types of 
spinons. By inserting (12) into (9), and using the formula 
a = 2 J2 P =± [Ac + Ap] , we obtain critical exponents for 
the oscillating parts, 

a? - 2(&r + £rc) , = 2(£ C r + e cC ) ■ (14) 

The numerical results for af and aj are plotted as a 
function of the magnetic field in Fig. 4. In weak mag- 
netic fields, represents how dominant the antiferro- 
magnetic fluctuations are, because pp ~ n/2. Since the 
value of a, is close to unity in weak fields, we can see 
that the antiferromagnetic correlation is indeed domi- 
nant there, but is gradually suppressed with the increase 
of the magnetic filed. For H > He, spinons of the C- 
sector become massive, and remaining massless spinons 
in the R-sector solely determine the critical behavior of 
the system. Accordingly, the third term of the right- 
hand side in eq.(13) vanishes. Also, the critical exponent 
af^ exhibits a discontinuity at H = He- This decrease in 
the critical exponent may imply that the correspond- 
ing spin correlation may be enhanced to some extent in 
magnetic fields. Similar phenomena have been indeed 
observed in a different context where spin correlations 
are enhanced by decreasing the multiplicity of bandsE3 . 
As seen from Fig. 4, however, the corresponding exponent 
is larger than unity, so this fluctuation may not become 
dominant. We note here that a logarithmic correction 
should enter the correlation functions at H = since the 
symmetry of the system is enhanced from U (1) to SU (2) 
at H = 0, although it does not appear in the present case 
with finite magnetic fields. _ 

Here, following the arguments of SachdevEJ) , we would 
like to briefly discuss the temperature dependence of nu- 
clear magnetic relaxation rates which may be an experi- 
mental probe to study low-energy properties of the sys- 
tem. Since the low-tcmpcrature behavior of the model is 
mainly determined by spinons mentioned above, we can 
obtain the leading temperature dependence of the nu- 
clear spin-lattice relaxation rates 1/Tj and the nuclear 
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C 3 T , (15) 
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Combining these expressions with the numerical results 
for the critical exponents shown in Fig. 4, we see that 
\jT\ behaves like ~ l/T 1 ^"^ at low temperatures for 
H < He- However for H > He this term vanishes and 
the low temperature behavior of 1/T% is determined by 
the term ~ 1/T 1_Q = . The temperature dependence of 
l/T^G f° r H < He is also governed by spinons in the 
C-sector and given by 1/T 3 / 2 ~ Q ? . For H > H c l/T 2G 
behaves like ~ 1/T 3 / 2 ~ Q ?. 

We have been concerned with the rather special model 
with alternating spins in this paper. Before concluding 
the paper, we wish to briefly mention what would be ex- 
pected if the system were away from the solvable point 
with coupling constants being changed. Unfortunately 
there have not been systematic attempts to this problem 
so far, so that we summarize here two possible disor- 
dered ground states predicted by the present analysis, 
(i) A possibility is that one of the c = 1 fields becomes 
massive in the presence of some relevant perturbations. 
In this case, we can see from the present results that 
low-field properties are classified by c = 1 CFT, and 
then the magnetization begins linearly. What is char- 
acteristic for this case is that a cusp structure in Fig.l 
may not be observed, (ii) Another interesting possibil- 
ity may be that both of two spinons become massive, 
and hence the system shows properties like a Haldane- 
gap system, for which there is no magnetization in low 
fields. These states, as well as the present massless state, 
should compete with the magnetically ordered state. In 
order to fully understand these problems, an extensive 
study with the aid of bosonization methods and numer- 
ical methods would be quite helpful. This issue is now 
under consideration. 

In summary we have investigated magnetic and criti- 
cal properties of an integrable alternating spin chain in 
external magnetic fields. The critical properties in mag- 
netic fields are classified in terms of CFT. Although the 
present analysis is based on a rather special model, it 
would provide a plausible guideline to systematically ex- 
plore characteristic properties of this new class of quan- 
tum spin chains. 



[1] G. T. Yee, J. M. Manriquez, D. A Dixon, R. S. McLean, D. M. 

Groski, R. B. Flippen, K. S. Narayan, A. J. Epstein, and J. 

S. Miller: Adv. Mater. 3 (1991) 309; Inorg. Chem. 22 (1983) 

2324; ibid. 26 (1987) 138. 
[2] H. J. de Vega and F. Woynarovich: J. Phys. A 25 (1992) 

4499. 

[3] H. J. de Vega, L. Mezincescu, and R. I. Nepomechie: Phys. 

Rev. B 49 (1994) 13223; Int. J. Mod. Phys. B 8 (1994) 3473). 
[4] S. R. Adams and M. J. Martins: J. Phys. A 26 (1993) L529. 
[5] M. J. Martins: J. Phys. A 26 (1993) 7301. 
[6] The cut-off parameters A^ are determined so as to minimize 

the free energy. 

[7] We note that a similar cusp structure has been found for the 



4 



M. Fujii, S. Fujimoto and N. Kawakami 



related spin model in which different spins are randomly dis- 
tributed. Sec P. Schlottmann, Phys. Rev. B 49 (1994) 9202. 
[8] H. J. de Vega and F. Woynarovich: Nucl. Phys. B 251 (1985) 

439 (1985); F. Woynarovich, J. Phys. A 22 (1989) 4243. 
[9] J.L. Cardy: Nucl. Phys. B 270[FS16] (1986) 186; H. W. Brote, 
J. L. Cardy and M. P. Nightingale: Phys. Rev. Lett. 56 (1986) 
742; I. Affleck, Phys. Rev. Lett. 56 (1986) 746. 

[10] A.G. Izergin, V.E. Korcpin and N.Y. Reshetikhin: J. Phys. 
A22 (1989) 2615. 

[11] H.Frahm and V.E. Korcpin: Phys. Rev. B. 42, (1990) 10553; 
N. Kawakami and S.K. Yang: Phys. Rev. Lett. 65 (1990) 2309. 

[12] T. Itakura and N. Kawakami: J. Phys. Soc. Jpn. 64 (1995) 
2321. 

[13] S. Sachdcv: Phys. Rev. B 50 (1994) 13006. 



Fig. 1. Magnetization as a function of the magnetic field. 



Fig. 2. Specific- heat coefficient as a function of the magnetic 
field. 



Fig. 3. Velocities of spin excitations as a function of the magnetic 
field. 



Fig. 4. Critical exponents af- and ajp for oscillating sectors in 
the spin correlator as a function of the magnetic field. 



